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Abstract: Agriculture needs to produce more food to feed the growing population in the 21“ century. It 
makes the reference crop water requirement (WREQ) a major challenge especially in regions with limited 
water and high water demand. Iran, with large climatic variability, is experiencing a serious water crisis due 
to limited water resources and inefficient agriculture. In order to overcome the issue of uneven 
distribution of weather stations, gridded Climatic Research Unit (CRU) data was applied to analyze the 
changes in potential evapotranspiration (PET), effective precipitation (EFFPRE) and WREQ. Validation 
of data using и situ observation showed an acceptable performance of CRU in Iran. Changes in PET, 
EFFPRE and WREQ were analyzed in two 30-а periods 1957—1986 and 1987—2016. Comparing two 
periods showed an increase in PET and WREQ in regions extended from the southwest to northeast and 
а decrease in the southeast, more significant in summer and spring. However, EFFPRE decreased in the 
southeast, northeast, and northwest, especially in winter and spring. Analysis of annual trends revealed an 
upward trend in PET (14.32 mm/decade) and WREQ (25.50 mm/decade), but a downward trend in 
EFFPRE (-11.8 mm/decade) over the second period. Changes in PET, EFFPRE and WREQ in winter 
have the impact on the annual trend. Among climate variables, WREQ showed a significant correlation 
(70.59) with minimum temperature. The increase in WREQ and decrease in EFFPRE would exacerbate 
the agricultural water crisis in Iran. With all changes in PET and WREQ, immediate actions are needed to 
address the challenges in agriculture and adapt to the changing climate. 
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1 Introduction 


According to the Food and Agricultural Organization (FAO), more than 3096 of the 
world's population depends on agriculture for livelihood in 2011. Water resources and land are 
necessary requirements for agriculture. Considering high demand for irrigation and water 
scarcity, water management and water conservation become very important. ш 1977, FAO 
developed a guideline to calculate the reference crop water requirement (WREQ) that was 
defined as the required water 
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to meet the consumed water by plants and compensate evapotranspiration (ET). The WREQ 15 
highly dependent on ET and effective precipitation (EFFPRE). 

ET, i.e., the sum of evaporation from canopy interception, bare soil evaporation and 
transpiration from vegetation, is a major component of water, energy and carbon cycle balance 
(Seneviratne et al., 2010). Global land ET returns about 60% of land precipitation into the 
atmosphere (Oki and Kanae, 2006). Besides, a complete space-time ET knowledge is critical for 
estimating the agricultural water requirement, as the world's largest water user (Beyazgül et al., 
2000; Hobbins et al., 2004; Xu et al., 2006; Burn and Hesch, 2007; Bouwer et al., 2008; Zhang et 
al, 2010; Liang et al, 2010; Espadafor et al, 2011; Sadeghi et al., 2015). Potential 
evapotranspiration (PET) is simply defined as the rate that ET occurs while the surface is 
well-supplied with water. Penman-Monteith (Allen et al., 1998) is a commonly used PET 
estimator with two variants equations for open water surface and idealized reference crop. Both 
equations calculate the PET based on the net radiation at the surface (R,), maximum and 
minimum temperature, wind speed and relative humidity (Allen et al., 1998). Any changes in 
individual climate variables would change the PET (Vicente-Serrano et al., 2014). Wu et al. 
(2019) analyzed the spatiotemporal variation in irrigation water requirement for wheat using 
effective precipitation and crop evapotranspiration. They calculated the irrigation water 
requirement which is different from the crop water requirement, which highly depends on the 
availability of meteorological data. 

Intergovernmental Panel on Climate Change (IPCC) reported changes in surface solar 
radiation, temperature, humidity and precipitation (Hartmann et al., 2013). Observational analysis 
using diverse measures suggested a drying of Earth's land in recent decades (Feng and Fu, 2013; 
McCabe and Wolock, 2015; Sheffield et al., 2012), consistent with precipitation and PET increase 
(Feng and Fu, 2013; Fu and Feng, 2014) and inconsistent with pan evaporation decrease 
(Roderick et al., 2015). A ubiquitous rise of estimated PET, driven by global warming (Scheff and 
Frierson, 2014), underlies the drying trends (Dai, 2011; Fu and Feng, 2014; РгауаПе, 2016; 
РгауаПе and Bandoc, 2019). 

PET is expected to increase directly as a result of surface heating (Trenberth et al., 2013), 
previous studies have analyzed the changes of ET globally and regionally. Globally an increase of 
PET anomalies has been reported from 1982 to 2013 (Zhang et al., 2015). The global trend of 
PET was significantly positive in the northeast and northwest of the United States, southern 
Africa, western Asia and coastal regions of Australia (Jung et al., 2010). 

In addition to global studies (Zhang et al., 2015; Zhang et al., 2016), the regional changes of 
evaporation have been analyzed in many countries such as United States (Walter et al., 2004; 
Irmak et al., 2012), China (Thomas, 2000; Wang et al., 2007; Zhang et al., 2007; Song et al., 
2009), India (Chattopadhyay and Hulme, 1997; Bandyopadhyay et al., 2009), Romania (Pravalie 
and Bandoc, 2015; Pravalie et al., 2019, 2020) and Iran (Tabari et al., 2011; Shadmani et al., 
2012; Tabari et al., 2012; Talaee et al., 2014). 

Previous studies associated the increase of PET to the relative humidity increase and solar 
radiation decrease (Chattopadhyay and Hulme 1997), an increase of temperature and solar 
radiation, and a decrease of relative humidity (Espadafor et al., 2011, Madhu et al., 2015). 
However, some studies found a decrease in PET in China that was associated with a significant 
negative trend of solar radiation and wind speed (Xu et al., 2006). 

A reliable analysis of PET requires accessible meteorological data which are not easily 
accessible especially in developing countries (Sherwood and Huber, 2010; Eyshi Rezaei et al., 
2015; Lobell et al., 2015; Zheng et al., 2012, 2015; Webber et al., 2017, 2018). Dinpashoh et al. 
(2011) obtained data for 16 stations in Iran and analyzed the changes of PET in 1965—2005. The 
changes did not show any specific pattern which can be explained by the low number of data. 
Wind speed was the most important parameter influencing PET (Dinpashoh et al., 2011). 
Analyzing PET in the western regions of Iran showed a positive trend in the majority of stations 
because of the increase of temperature during the study period (Tabari et al., 2011). The 41-a 
(1965—2005) data were obtained from 11 stations to study the changes of PET in arid regions of 
Iran (Shadmani et al, 2012). Limited access to long-term, quality-controlled and 
spatially-distributed data encourage applying the gridded global/regional datasets as promising 
alternatives to in situ observations. 

In addition to PET, EFFPRE 15 another important parameter to calculate the WREQ. In 
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definition, EFFPRE is the amount of rainfall that is actually added to the soil and can be used by 
plants. The EFFPRE does not include the parts of rainfall that flow over soil (runoff) and 
percolates below the root zone (deep percolation) of plants (Dastane, 1974). Parameters including 
climate, depth of the root zone, soil structure, initial soil moisture content, topography and soil 
texture influence the EFFPRE (Brouwer and Heibloem, 1986). EFFPRE is a crucial factor when 
designing and operating irrigation systems and can be calculated based on different methods 
(Dastane, 1974; Brouwer and Heibloem, 1986; Patwardhan et al., 1990) depending on the 
simplification of the hydrologic cycle. EFFPRE has been calculated for diverse applications using 
different data (Karamouz et al., 2004; Sanaee-Jahromi et al., 2000; Kashani et al., 2016). 

Agriculture plays a significant role in Iran's economy. Agriculture is the major source of food 
and contributes an average of 20% to the gross domestic product. It provides one-third of the 
country's jobs and 90% of agricultural raw materials and accounts for half of the total non-oil 
exports. Considering the limited source of water in most parts of Iran, most regions have always 
suffered from an inefficient agricultural system that relies highly on irrigation (Nattagh, 1986). 
How PET, EFFPRE and МКЕО are changing over time in a changing climate, would seriously 
impact Iran's economy. Here, the interest was to calculate the EFFPRE with complete spatial 
coverage and analyze the changes of EFFPRE over time. Considering the limited data access in 
Iran, 30-а meteorological parameters were obtained from the Climatic Research Unit gridded 
Time Series (CRU TS) dataset to investigate spatial changes of ET, EFFPRE and WREQ over 
time. 


2 Materials and methods 
2.1 Study area 


Iran, located in southwest Asia, is a mostly mountainous country with a topography ranging from 
—28 m (on the southern coast of the Caspian Sea) to 5671 m (in Mount Damavand). With an area 
of 1.6x10° km’, Iran is the second-largest and the second most populated country in the Middle 
East. The climate of the country 15 highly influenced by two major mountain chains in the north 
(Alborz) and West (Zagros), the Caspian Sea in the north, and the Oman Sea and Persian Gulf in 
the south. Iran has a diverse climate from arid (6596) and semi-arid (2096) to humid and 
semi-humid (15%) climate (Madani, 2014; Khalili et al., 2016). January and July are the coldest 
and hottest months of a year with temperatures ranging —6?C—21?C and 19?C—39?C, respectively 
(Madani, 2014). Diverse weather systems are discovered in Iran. In the northwest, cold winters 
with heavy snowfall are expected in December and January and the weather is relatively mild 
during spring and autumn. However, in the south and southeast, moderate winters and hot-humid 
summers are expected (Saboohi et al., 2012). Precipitation 1s highly variable in space and time 
(Raziei et al., 2012, 2014). Annual precipitation ranges from 50 mm in deserts of central Iran to 
1000 mm in the southern coasts of the Caspian Sea (Madani, 2014, Валет et al., 2014). About 
70% of the precipitation falls between November and March (Raziei et al., 2005) over only 2596 
of the region (Madani, 2014). The annual mean minimum temperature (TMN) and maximum 
temperature (TMX), and annual total precipitation in two 30-year periods of 1957—1986 and 
1987—2016 are showed in Figure 1. 

Almost 3.7x10? km? of the countries surface area are arable lands (Abbaspour et al., 2009) 
including 1.9x10? km? devoted to field crop and horticulture (Keshavarz et al., 2005). More than 
half of the lands (1.0x10° кт”) are rain-fed and the rest are irrigated applying traditional or 
modern techniques (Keshavarz et al., 2005). Only 5% of irrigated lands in Iran are under 
pressured irrigation and the average irrigation efficiency is 35% (Madani et al., 2014). 


2.2 Data 


The CRU TS 4.01 dataset (Harris et al., 2014) has been developed by the UK's Natural 
Environment Research Council (NERC) and the US Department of Energy to build a global 
climatic database. The well-tested, popular and publicly available gridded data contains 
information for 10 different climate variables, interpolated from 5000 stations ш 1901—2016. 
Spatially, data covers the entire globe between 6095 and 80°N (excluding Antarctica) 
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Fig. 1 Annual mean minimum temperature (ТММ), maximum temperature (TMX), and total precipitation (PRE) 
in two 30-year periods of 1957—1986 and 1987—2016. The panel on the right shows the difference between the 
two periods. Grid-cells with non-significant changes (90% confidence level) are shown in grey. 


at 0.5? resolution. Harris et al. (2014) compared the CRU TS dataset with other similar gridded 
datasets and showed a significant agreement between datasets and sparse observations. Previous 
studies validated the CRU dataset especially for trend analysis (Harris et al., 2014; Shi et al., 
2017). Mean monthly temperature and PET, and the monthly total precipitation was obtained for 
675 grid-cells for 1957—2016 in Iran (Fig. 2). PET has been calculated from Penman-Monteith 
equations using gridded temperature, vapor pressure and cloud cover data (Harris et al., 2014). 


2.3 Definitions 
EFFPRE in this study is estimated based on the average of the United States Department of 
Agriculture (USDA, Eq. 1) and FAO (Eq. 2) methods which have been validated previously for 


Iran (Rahimi et al., 2013). Both methods simplified the water balance into just using precipitation 
when calculating the EFFPRE. 


Px(125.0.2 x P) 
EFFPRE = 125 
1254 0.1P P > 250 mm 


P «250 mm 
А (1) 
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0.6P, -10 P, € 70 mm 
EFFPRE - , (2) 
0.8P, -24 P >70 шт 


where P and Р„ are the annual and monthly precipitation (mm), respectively. Once the EFFPRE is 
calculated, the WREQ would be calculated based on the difference between PET and EFFPRE. 
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Fig. 2 Selected grid-cells from the СКО TS (Climatic Research Unit gridded Time Series) dataset and 14 
stations used to validate the dataset 


2.4 Statistical test 


Two different sets of statistical methods were used in this study. First, linear parametric Pearson 
correlation (r) and root mean square error (RMSE) were used to validate the CRU TS data in Iran. 
These methods have been widely used to validate diverse sets of data worldwide (Entekhabi et al., 
2010; Tavakol et al., 2019). Second, the ranked-based, non-parametric and distribution-free 
Mann-Kendall (MK) test (Mann, 1945; Kendall, 1975) has been used to analyze the temporal 
monotonic trend of variables. Besides, the trend Free Pre-whitening (TFPW) procedure was 
applied to remove the effect of serial correlation (Yue et al., 2002). The R package 'zyp' was used 
to carry out the test (Bronaugh and Werner, 2015). The confidence level for the tests was 90%. 


3 Results and discussion 


3.1 Validation of CRU TS data using in situ observations 


Considering the importance of temperature and precipitation on both PET and EFFPRE, the 
accuracy of CRU TS data was evaluated against observed data obtained from 14 selected stations 
for 1951-2016. For monthly temperature (Fig. 3), the average А? and RMSE between CRU TS 
and observed data are 0.97 ?C/decade and 2.17?C/decade, respectively. Except for Anzali station, 
other stations showed acceptable results. For monthly precipitation (Fig. 4), the average №? and 
RMSE between CRU TS and observed data are 0.85 and 10.71 mm/decade, respectively, 
excluding the Anzali station. The results of this validation confirmed the findings of the previous 
study on the acceptable performance of CRU TS data in Iran (Miri et al., 2016). 

Due to the attempt of this study to analyze the trend over time, the accuracy of the trend 
estimated by CRU TS data was validated before and after 1987, separately (Fig. 5). The 
performance of CRU TS precipitation data in estimating the trend was improved for the period 
before 1987 (RMSE, 3.23 mm/decade) to the period after 1987 (RMSE, 2.27 mm/decade). For 
monthly temperature data, the same improvement was found in estimating the trend before 
(RMSE, 0.56?C/decade) and after 1987 (RMSE, 0.18?C/decade). CRU TS data showed 
acceptable results calculating the trend for both precipitation and temperature, except in the 
Caspian coast in the north. Overall, the CRU TS data showed a better performance in the recent 
30-a period which can be explained by the increase in the number of observed stations as well as 
the data quality (Harris et al., 2014). 
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Fig. 3 Relationship between CRU TS data and observed monthly temperature data in 14 selected stations and all 


stations 
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Fig. 4 Relationship between CRU TS data and observed monthly precipitation data in 14 selected stations and 
all stations. Except for Anzali, most stations showed a good agreement between datasets. 
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Fig. 5 Relationship between trends calculated from CRU TS and observed monthly precipitation (left) and 
temperature (right) data before and after 1987. Acc, accuracy; Decr, decrease; Incr, increase. 


3.2 Changes in long-term averages 


321 PET 


Figure 6 shows the mean annual and seasonal PET in two 30-a periods before and after 1987. 
North and northeast regions showed no significant difference for annual PET between two periods 
(the first row in Fig. 6). In the southeast, a decrease was discovered in 1987—2016 compared to 
1957—2016. The decrease was significant in all seasons. However, a significant increase in PET 
was found in the east, southeast, southwest and central regions of Iran. The highest increase in 
PET occurred in spring and then summer (Fig. 6). However, the changes in PET were not 
significant in these regions in winter. Monthly analysis of PET value showed the highest values of 
PET in June, July and August. Spatially, the highest PET was recorded in the east and west during 
these months. 


3.2.2 EFFPRE 


The mean annual and seasonal EFFPRE in two 30-a periods are presented in Figure 7. Overall, 
the arid and semi-arid regions of Iran are surrounded by major mountain chains that extend from 
central to the southwest. These regions had the lowest values of EFFPRE. However, Zagros 
Mountains extending from west and northwest to southwest showed the highest EFFPRE. 
Annually, most regions showed non-significant change and the regions with significant decrease 
were in the northwest, northeast, and southeast (Fig. 7). Seasonally, the highest and lowest 
EFFPRE were found in winter and summer, respectively. The difference between the two periods 
is non-significant in summer and autumn. In spring and winter, the decrease occurred in the 
southeast, northwest and northeast (Fig. 7). Agriculture is a very important sector of the economy 
in these three regions. Normally, Iran receives the least precipitation in summer and any change 
of EFFPRE in summer has no significant influence on agriculture. However, the decrease of 
EFFPRE in winter and spring in the northeast, northwest and southeast has a crucial negative 
influence on the crop yields especially for rainfed crops and autumn-planted grains. In the 
northwest, the grains such as wheat and barley are mainly planted in November for early summer 
harvest in June after receiving winter precipitation. Legumes including lentils and peas are also 
planted in March to be harvested in summer. In the southeast, most plants are planted in March to 
benefit the spring and winter precipitation. 

3.2.3 WREQ 

WREQ has been calculated based on the difference between PET and EFFPRE. The spatial mean 
annual and seasonal WREQ almost follow the same pattern of PET. Mean 30-a values of WREQ 
increase from the northwest to southeast for both periods. Annually, an increase was found in the 
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Fig. 6 Annual (top panel) and seasonal (bottom panel) averages of total potential evapotranspiration (PET) in 
two 30-year periods of 1957-1986 and 1987—2016. The panel on the right shows the difference between the two 
periods. Grid-cells with non-significant changes (9096 confidence level) are shown in grey. 
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Fig.7 Annual (top panel) and seasonal (bottom panel) averages of total effective precipitation (EFFPRE) in two 
30-year period of 1957—1986 and 1987-2016. The panel on the right shows the difference between the two 
periods. Grid-cells with non-significant changes (9096 confidence level) are shown in grey. 
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south, east and northeast during 1987—2016 compared to 1957—1987 (Fig. 8). Seasonally, the 
highest increase in WREQ happened in spring and then summer. However, a decrease in WREQ 
was found in the southeast in summer. The increase of WREQ in summer and spring can be 
explained by the increase of PET in the same regions. In northeast, EFFPRE played a crucial role 
in the increase of WREQ, especially in spring. The mean seasonal WREQ is higher in summer in 
1987—2016 compared to 1957—1987 in the southeast and northwest due to PET drop. Between 
April and September, most regions showed very high WREQ, making irrigation necessary. 

The increase of WREQ in arid and semi-arid regions in the center, south, east, and northeast 
would cause major losses in groundwater storage because of over-pumping in spring and summer. 
Considering the increase of WREQ in the second 30-a period, discovered mainly in June, July 
and August, there would be a need for better management with an efficient irrigation system. No 
significant increase of WREQ in the second 30-year period in autumn and winter encourages 
more planting during autumn and winter. 


3.3 Trends 


3.3.1 РЕТ 


The MK test was used to analyze the trend of annual РЕТ in Iran. In the first period, an upward 
trend was discovered in the south and southeast (Fig. 9). However, there was a downward trend in 
the central northeast. In the second period, the trend was mostly positive (6196 significant and 
38% non-significant) in most regions except in areas extended from southwest to northeast. The 
slope of the trend showed a decrease in the south, but an increase in west, northwest, central 
northwest and east. Hosseinzadeh et al. (2013) showed an increase of PET in seven stations out of 
4] stations during 1966—2005. The stations with significant upward trend were mainly located in 
the northeast and northwest. Tabari et al. (2012) discovered an increase in monthly values of PET 
from January to July and a decreased from August to December for 21 weather stations across the 
country. 

Tabari et al. (2012) estimated an average trend of 0.44 mm/a in annual PET in the last four 
decades. The value was 70% smaller than the value estimated in this study. It can be explained 
with the different nature of data and diverse study periods. The location of stations in the previous 
study is not well-distributed in the country and could not cover the northern regions where a 
significant increase was detected in 1987—2016. 

Seasonally, winter showed a considerable difference between the two periods. In 1957—1986, 
there was a significant negative trend in the northern half of Iran. In 1987—2016, however, the 
trend was mostly positive in the majority of regions (89%) except the northeast (Fig. 9). Talaee et 
al. (2014) suggested the most significant trend of PET in winter (1966—2005) considering the 
magnitude of trend and the frequency of stations with a significant positive trend. The average 
slope of the positive trend was between 3.17 and 14.32 mm/decade over the first and second 30-a 
periods, respectively. The significant trend in spring and summer was mainly positive and mostly 
discovered in the south, east and southeast. In the autumn, the trend was negative in the west and 
northwest in the first period. However, in the second period, a large area of northeast (24%) 
showed a negative trend as well (Fig. 9). Monthly trend analysis of PET showed a significant 
upward trend in the majority of the area in March, June and January. However, the majority of the 
area showed a downward trend in November. For other months, the trend was either 
non-significant or positive. 

3.3.2 EFFPRE 

Annually, the first period (1957—1986) showed an increase (10.6 mm/decade) in the EFFPRE in 
the Zagros Mountains toward the north. However, the second period showed a decrease in the 
regions towards west and southwest of the Zagros Mountains, east and central northeast Iran (Fig. 
10). No decrease has been recorded in EFFPRE in the first period. In 1987—2016, most areas 
(93%) showed a significant (42%) and non-significant (5196) negative trend with an average of 
11.8 mm/a (Fig. 10). The decrease would be explained by the decrease of EFFPRE in winter. 

Seasonally, the trend was non-significant for most of the seasons in both periods except for 
winter and autumn. In winter, the trend was increasing and decreasing in the first and second 
periods, respectively. The first period showed an upward trend in autumn. Monthly, the trend was 
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Fig. 8 Annual (top panel) and seasonal (bottom panel) averages of water requirement (WREQ) in two 30-year 
period of 1957-1986 and 1987-2016. The panel on the right shows the difference between the two periods. 
Grid-cells with non-significant changes (9096 confidence level) are shown in grey. 
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Fig. 9 Temporal trends in annual (top panel) and seasonal (bottom panel) total potential evapotranspiration 
(PET) in 1957-1986 and 1987—2016. Grey areas show non-significant trends. Panels on the right represent the 
difference between the slope of trends in two 30-year periods. 
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significantly negative in a few grid-cells in January, February, March and December. About 60% 
of the region showed an upward trend in EFFPRE in November. 

Gu and Adler (2015) reported a decrease in precipitation in 1979—2012 globally. In Iran, 
analysis of the NASA GISS Climate Model showed a negative trend (—2.6 to —6.5 mm/a) in 
1951—2009 (Raziei et al., 2014b). 


3.3.5 WREQ 


Annually, WREQ was significantly decreased in the west and central northwest and was 
significantly increased in the south and southeast in the first period (Fig. 11). That can be 
explained by the drop of PET and the rise of EFFPRE at the same time. In the second period, 
most regions experienced a significant (5896) and non-significant (40%) increase in the WREQ 
with an average of 25.5 mm/decade. Winter represented the highest difference between the first 
and second periods. The WREQ was decreasing in the first period and increased in the second 
period. The trend was mostly upward in spring and summer. However, the trend was mostly 
decreasing in autumn. The temporal trend of WREQ follows the same pattern of PET in winter, 
summer, and spring and the same pattern of EFFPRE in autumn. 

Spring and summer are the most important agricultural seasons in Iran. Irrigation starts usually 
in April and continues to September. Results of this study showed an increase in the WREQ in 
spring and summer in 80% and 93% of areas, respectively. The WREQ rise is more important in 
the south and southeast with a significant water crisis. These regions are covered with large-size 
gardens of citrus, pistachio and dates. People heavily rely on agriculture in these regions, 
especially in summer. The monthly analysis showed a significant decrease in November and a 
significant increase in other months. 


3.4 Correlation analysis 


The spatial correlation between defined variables was tested using the Pearson correlation test. 
The correlation coefficient changed from —0.6 between annual EFFPRE and МКЕО to 0.7 
between PET and WREQ. As expected the mean annual WREQ increase when PET increases. 
However, the annual WREQ decrease when the EFFPRE increases. Seasonally, the greatest 
association between EFFPRE and WREQ was discovered in spring (r= –0.92), winter (r= —0.82) 
and summer (r= —0.73). In autumn, the correlation between EFFPRE and WREQ was less strong 
("= —0.44) compared to other seasons. However, the association between PET and WREQ was 
greater in autumn (7—0.66), spring (7=0.64) and winter (7=0.58). Summer showed a lower 
correlation (7=0.43) between PET and WREQ. Results showed the more significant role of 
EFFPRE on WREQ in spring, winter and summer. In autumn, PET plays a more important role in 
WREQ than EFFPRE. 

Climate variables including mean temperature, maximum temperature, minimum temperature 
and precipitation influence PET, EFFPRE and WREQ. Correlation analysis showed a positive 
correlation between PET and minimum temperature (7=0.52), mean temperature (7=0.63) and 
maximum temperature (7=0.54). However, a negative relationship was discovered between PET 
and precipitation (7= —0.53). The same pattern of association was discovered between WREQ and 
minimum temperature (7=0.59), mean temperature (7=0.46), maximum temperature (7=0.15) and 
precipitation (7= —0.29). Except for the minimum temperature, the correlation between climate 
variables and EFFPRE is lower than the correlation between climate variables and PET. It can be 
explained with the significant correlation (r= —0.5) between minimum temperature and 
precipitation. The decrease in precipitation in Iran can explain the increase in minimum 
temperature. 


4 Conclusions 


Determination of PET and WREQ is necessary for agricultural water management, especially in 
arid and semi-arid regions. Changes in climate variables such as temperature, precipitation, wind 
speed and relative humidity are supposed to change the WREQ. Long-term weather data are 
required to calculate WREQ in a specific geographical region. Due to limited observation data in 
Iran, the CRU TS gridded data has been used to investigate the spatiotemporal changes of PET, 
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Fig. 11 Temporal trends in annual (top panel) and seasonal (bottom panel) water requirement (WREQ) in 
1957-1986and 1987—2016. Grey areas show non-significant trends. Panels on the right represent the difference 
between the slope of trends in two 30-year periods. 
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EFFPRE and WREQ. Validation analysis of data showed the acceptable performance of data in all 
stations except for the Anzali station located in the north. 

Results showed a decrease of EFFPRE with an average rate of 11.8 mm/decade in the second 
30-a period (1987—2016). This confirmed the results of previous studies on the precipitation drop 
in Iran. An increasing trend in PET (14.32 mm/decade) and WREQ (25.5 mm/decade) was also 
discovered in the second 30-a period. The increase was more significant in winter. 

In the second period (1987—2016), the increasing trend of PET in winter and especially in 
March showed a shift in the growing season. The earlier beginning of the growing season is very 
important for autumn season farming crops including wheat, rapeseed and maize. In addition, a 
decreasing trend of PET especially in November showed the early beginning of the freezing 
season in all regions. The shift of growing season may suggest a shift in the cultivation of autumn 
crop farming toward October before the beginning of the freezing season. For EFFPRE, a shift 
was seen in the rainfall season. The rainfall season began earlier in the second period compared to 
the first period. Results also showed an increase in EFFPRE in November. Shifting the irrigation 
period to the months with the highest EFFPRE can help water management especially in regions 
facing a water crisis. Results showed an increase in EFFPRE since March. So, shifting the 
irrigation to March from April gives enough time for plants to grow before June when the weather 
is too hot and irrigation consumes much water. Considering the influence of EFFPRE drop in the 
second period during winter, encourage more management in agriculture. 

Agriculture in Iran is completely dependent on irrigation. In Iran, about 90% of water resources 
are used for irrigation. Results showed either a significant or non-significant increase in WREQ in 
the majority (98%) of areas. The increase of WREQ suggests a probable pressure on groundwater 
resources for Iran with limited surface water supplies. The significant decrease of groundwater in 
Iran confirmed the impact of irrigation. Modifying agricultural patterns, using the plants resistant 
to drought, and more efficient irrigation systems are some suggested strategies to cope with the 
water crisis in Iran. There is also a need to update the national WREQ for Iran considering the 
changing climate. 
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